The Ground State and Excitations of Spin Chain with Orbital Degeneracy 
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The one dimensional Heisenberg model in the presence of orbital degeneracy is studied at the 
SU{4) symmetric point by means of Bethe ansatz. Following Sutherland's previous work on an 
equivalent model, we discuss the ground state and the low-lying excitations more extensively in 
connection to the spin systems with orbital degeneracy. We show explicitly that the ground state 
is a 517(4) singlet. We study the degeneracies of the elementary excitations and the spectra of the 
generalized magnons consisting of these excitations. We also discuss the complex 2-strings in the 
context of the Bethe ansatz solutions. 
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I. INTRODUCTION 



There has been much interest recently in spin Hamil- 
tonians with orbital degeneracy. The orbital degree of 
freedom may be relevant to many transitional metal ox- 
ides ]l[-pd]|. Examples of such systems in 1-dimension 
include quasi-one-dimensional tetrahis-dimethylamino- 
ethylene_(TDAE)-C6o P^l, and artificial quantum dot 
arrays [Q. Recently, we discussed these systems jlj] 
within the framework of a S'C/(4) theory. A quantum 
disordered ground state in 2-dimension was proposed to 
be relevant to the experimentally observed unusual mag- 
netic properties in LiNi02- There have also been nu- 
merical studies of the 1-dimcnsional models for these sys- 
tems ffsl 



16|. 



In the present paper, we use the Bethe ansatz method 
to study the S'C/(4) symmetric Heisenberg spin chain with 
two-fold orbital degeneracy. This model is equivalent to 
the model studied by Pokrovskii and Uimin p7[ , and to 
one of a class of models that has been solved by Suther- 
land pq ]. Expanding on Sutherland's work, we study the 
ground state and low-lying excitations more extensively 
by considering holes and 2-strings in the thermodynamics 
limit, and in connection to the spin systems with orbital 
degeneracy. We show explicitly that the ground state 
is a SU{4:) singlet, consistent with the generahzed Lieb- 
Mattis theorem of Affleck and Lieb |g^ . We discuss the 
degeneracies of the elementary excitations and the spec- 
tra of the generalized magnons resulting from such such 
excitations. We also discuss the complex 2-strings in the 
context of the Bethe ansatz solutions. The paper is or- 
ganized as follows. In section II, we introduce the model 
and discuss its symmetry. In section III, we present the 
Bethe ansatz solution following Sutherland's approach. 
We discuss the ground state in section IV, and the ele- 
mentary excitations and the generalized magnon modes 
in section V. A brief summary is given in section VI. 



II. SYMMETRY CONSIDERATION 

We consider the spin chain of N sites with two-fold 
orbital-degeneracy and with periodic boundary condition 



N 



n 






(2r,-f,+i + i)(25,.5,+i + i) 



(1) 



where Sj and Tj are the spin and orbital operators respec- 
tively on the j-th site, and are each generators of a SU{2) 
Lie algebra. Clearly, in addition to the permutation sym- 
metry (H is invariant under the transformation j —> j + 1 
for all j's), the system has an explicit SU{2) (g) SU(2) 
global symmetry and a bi-symmetry Sj <-^ Tj. Actually 
the Hamiltonian (|l|) is invariant under a global S'[/(4) 
transformation, which is generated by S" = ^2^^^' 

T" = EjTj": and F"^ = V . T^" ® S^ {a, (3 = 1,2,3). 
These operators satisfy the following commutation rela- 
tions 

[S°' , S^] = le'^'^^ S^ , [S'",T'']=0, 

[ya/.^ y/3^] ^ i^o^P-rgt^'^S'^ + ie'"'Pd°''^TP. (2) 

Such a symmetry was noticed by Wigner in the study of 
nuclei long time ago [^. Since the group SU{4) is of 
rank three, there are three conserved quantum numbers 
in general. It is useful to write the above commutation 
relations in terms of Chevalley basis, i.e., three genera- 
tors in the Cartan sub- algebra of SU{A) (precisely, the 
A3 Lie algebra) if„ = J2j=i Hnij) {n = 1, 2, 3) which 
can be diagonalised simultaneously, and the other twelve 
generators Ea = J2j=i^a{j) (a denotes root vectors). 
The local generators Hn{j) and Ea{j) are related to the 
spin and orbital operators by, 



H,{j) = (1 + 27;^)5|, E^, (j) = (- + T^)S+, 



H,ij)^iT^^S^ 



EcAj) 



nsj 



H,{j) = (1 - 2T^)S^, E^M = (- - T^)S+, 

-Eai+aa+QsO) = ^j '5'j 7 (3) 

where Qi, a2 and as denote the simple roots of ^3 Lie 
algebra |^l|. The generators corresponding to the nega- 
tive roots are given by E-a{j) = El^{j). One can verify 
that these operators satisfy the standard commutation 
relations of A3 Lie algebra, and Eq.(P can be rewritten 
as, 



where A denotes the set of roots of the Lie algebra. We 
denote the spin components by up (t) and down (|), and 
the orbital components by top and bottom. Then the 
four possible states on each site are 



|1> 
|2> 



:=|I >-|l/2,l/2>, 
:=|1 >=|-l/2,l/2>, 



|3> :=| I >=|l/2,-l/2>. 



14 > := 



i 



>=|-l/2,-l/2> 



The local lowering/raising operators £'±q„ (j) relate those 
four states on the jth site as follows. 



E^a„{j)\n >j -^ \n+l >j, 
Ea„ij)\n + 1 >j —^\n>j, n = 1,2,3. 



(4) 



hU,J + ^) = 5]5'""i^m(j)i?n(j + l)+E Eaij)E.^ij + l), 



qSA 



In terms of those operators, a general state can be written 
as 



M" 



M' 



M 



IV' >= Yl ^{x:y;z)Y[E^aAzk)'[[E^a2{y3)Y[E^ai{xt) 

{xi}{yj}{^k} k=l j=l i=l 



T T 



(5) 



where x := {xi,X2,- ■ ■ ,xm), V ■= ivi, y2,- ■■ ,yM'), 
z := {z-Y,Z2,--- ,ZM"):. 1 < xi < X2 < ■ ■ ■ < xm < N; 
xi < vi < y2 < ■ ■ ■ < Vm' < Xm, yi < zi < z2 < 
■ ■ ■ < ZM" < yM'\ and {zk} C {y,} C {xi}. We may de- 
fine the weights as the eigenvalues of the global operator 
Hn, indicated by {Hi, H2, H^). The eigenvalues of the lo- 
cal operators Hi{j), H2{j), H3{j) acting on the four local 



states 



T 



> 



i 



J, I _ >j,l I >3 and I I >j are (1,0,0), 

(-1,1,0), (0,-1,1) and (0,0,-1) respectively. We shall 
focus on the state with the highest weight. The other 
states in the same irreducible representation can then be 
obtained by using the corresponding lowering operators 
E-curi ■ In the present model, the irreducible represen- 
tation of the S'[/(4) group of a N-site system is labeled 
by 

(N + M' - 2M, M + M" - 2M', M' - 2M"). (6) 



III. THE BETHE ANSATZ SOLUTION 

The permutation and the S'C/(4) symmetries in the 
Hamiltonian enable us to seek the eigenstate of both the 



cyclic permutation operator and the generators of the 
Cartan subalgebra of A3. The invariance of the cyclic 
permutation imposes a periodic boundary condition on 
the wave function ip(x,y,z). The present model is solv- 
able [Q, and the Bethe ansatz equations for the spectra 



Xj + i/2 
J~i/2 



_ _TT ^j - M +1 TT fJ-p - Aj +1/2 

^^^■-\l~^^}^^^p-X,-^/2' 



1=1 



M 



-p-r ^-y - A; + t/2 _ _ -p-r ^^ - flp+l -p-r l^b - fJ'f + 1/2 

Ai ^^-\i- i/2 J-^i fi-^- fip-i ii i/fc - /i^, - i/2 ' 

M' .. M" . 



/3= 



\ i^c- fJ-fi - «/2 



b=l 



I'c- I'b -i 



(7) 



where j,l = 1, 2, • • • , M; /3, 7 = 1,2,---,M' and 6, c = 
1, 2, • • • , M" . These are secular equations for the spectra 
of 5/7(4) rapidities A, jj, and ly. The energy spectrum is 
given by 



M 



E ■ 



E 



J 



^(1/2)2 + A? 



2 ■ 



(8) 



The momentum defined by the translation of the system 
along the chain is given by 



1=1 ' 

M 

= ^[7r-2tan-i(2A0]. 
1=1 



(9) 



Note that P in eq. (g) is determined up to mod {2tt), 
and the inverse trigonometric function is defined in the 



main branch. We have included explicitly the tt term in 
eq.(|g|), which is usually neglected in the study of pure 
spin Heisenberg models [gO[. In the S't/(4) model, there 
are three types of elementary excitations as we will dis- 
cuss below, and it is convenient to include the tt term in 
P to study the magnon types of composite excitations. 
We define the momentum of the elementary excitations 
as the momentum relative to the ground state |20[| . By 
taking the logarithm of (R), a set of coupled transcenden- 
tal equations are obtained. 



M 



M' 



ei/2(A,) - ^Eei(A, - A,) - l^e,/2(M/3 - a,) = ^/„ 

1=1 (3 

M M' M" 

^61/2(^7 -A;) - ^Oiifi^ - fip) -^61/2(1^6-^7) = 27rJ^, 

1=1 P=l b=l 

M' M" 

y^ 6i/2(!^c - M/3) - X! ^l('^c - l^b) = 2ttKc, 

13=1 b=l 



(10) 



where 6^(2:) '■— 2tan~^(a;/p). The quantum number 
Ij is an integer or half-integer depending on whether 
N — M — M' is odd or even, and so is J^ (or Kc) de- 
pending on whether M - M' - M" (or M' - M") is odd 
or even. These properties arise from the logarithm func- 
tion. 

Replacing \j , /i-y and Vc in eq. (nOh by continuous vari- 
ables A, /i and v but keeping the summation still over 
the solution set of these roots {A;,/i^, I'fc}, we can con- 
sider the quantum numbers Jj, J^ and Kc as functions 
/(A), J(/i) and K{y) given by eq. (|^. Take /(A) as an 
example. When /(A) passes through one of the quantum 
numbers Ij , the corresponding A is equal to one of the 
roots \j. Similarly for J{p) or K{v). However, there may 
exist some integers or half-integers for which the corre- 
sponding A (/i or v) is not in the set of roots. We shall 
name such a state as a "hole". In the thermodynamics 
limit iV ^ 00, we may introduce the density of roots and 
the density of holes (indicated by a subscript h) , 



a{\) + an{\) = [l/N)dI{\)/dX, 

uJ{^l) + LJhi^i) - {l/N)dJ{^l)/dfl, 

T{v)+Th{v) = {l/N)dK{u)/dv. 

By replacing the summations by integrals, 

1 M B 

1=1 ''^^ 

^°° [3=1 ''^^' 

N^oo iV f— f J_B" 



6=1 



eq.(9) become the coupled integral equations. 



a{X) + <Jh{X) = i^i/2(A) - / dX'K,{X - A')a(A') + / dfi'Ky^iX - fi'Mfi'), 



B 

cuifi) + ujhifi) = I dA'i^i/2(A* - A')a(A') - 

-B 

r-B' 

t{v) + Th{y) = I d^'Ki/2{v - fi')uj{fj,') 

B' 



B' pB" 



dv'Ki{v~v')T{v'), 



(11) 



where isrp(a;) := tt ^p/{p'^+x'^), and B, B' , and B" inthe 
definite integrals should be determined self-consistently. 
In the absence of the complex roots, M/N — J_g a{X)dX, 

M'/N = J^g, uj{fi)dfi, and M"/N = J^'^,, T{v)dv. Once 
the density a is solved from eq.(|ll|), we have the z- 
components of the total spin and the total orbital 



S', 



tot 

N 



B' 



= -+ / Lo{p)dp— / (j{\)dX — 



T, 



tot 

N 



B' 
B' 



uj{p)dfj., 



T{v)dv, 



(12) 



the energy 

E = -2t:NJ 
and the momentum 



Ky^{\)a{\)d\ 



P = -N 



[2tan"i(2A)-7r]cr(A)dA. 



IV. THE GROUND STATE 

The ground state is described by the densities cro(A), 
Wo (a*), and tq{v) with no holes and by Bq = Bq = Bq -^ 
oo. This is true because all the states with holes will have 
higher energies. In this case, eq. ( [ll| ) can be solved. Let 



^o(A) - ^ 






1 

2^ 
1 

' 2^ 



aoiq)e-'''^dq, 
Cuo{q)e-'''^dq, 

) 

fo{q)e-''>''dq, 



then 



(Joiq) = sinh(3(7/2)/sinh(2q), 
uja{q) = sinh(q)/sinh(2q), 
To(g) = sinh((7/2)/sinh(2(7). 

Hence from eq. (^, the highest weight labeling the 
ground state is the null vector (0, 0, 0), and so the ground 
state is a SU{4) singlet. This agrees to the theorem of 
Affleck and Lieb [E2| , a generalization of the spin chain 
problem |2^. In that state, the total orbital, the total 
spin, and their products are all zero, i.e.. 



T, 



tot 

N 



_ 1 
~ 2 

CZ 
'-'tot 

N 



ujo{p)dp = - -i2}oi0) = 0, 
+ c^o(0)-ao(0)-fo(0) = 0, 

AT 



In deriving this result, eqs.(Q) have been used. The en- 
ergy and the momentum of the ground state are 



En — 



-7VJ(^ln2+J) 

mod 2n, for N/A — even, 
TT mod 27r, for N/A ~ odd. 



(13) 



Eq. ( [131) coincides with the result of Sutherland jTsj af- 
ter correcting for the trivial overall constant shift JN 
between the two models. 



V. LOW-LYING EXCITATIONS 

A. Spectra of elementary excitations 

The possible elementary excitation modes are obtained 
by the variation in the sequence of quantum numbers 
{Ij}, {J-y} or {Kc} from the ground state. We can 
assume B — B' = B" -^ oo for the low-lying exci- 
tations. The simple modes will be solved by placing 
holes in the rapidity-configurations. If we let a{X) = 
CTo(A) + CTi(A)/7V, uln) = wo(m) + wi(/.t)/A^ and T{iy) = 
To(i^) +Ti(iy)/N, then the excitation energy and momen- 
tum 



AE = -2ttJ / i^i/2(A)CTi(A)dA, 

J — OO 

/OO 
[2tan"i(2A)-7r]cri(A)dA, (14) 

-DO 

and AM = J ai{X)dX, AM' = /wi(^)d/x, and AM" = 
^ Ti{v)dv. After solving the integral equations ( pd| ) 
with ah{X) = 5{X - X)/N, w/i(^) = S{n - fl)/N or 
Th{i^) = S{i' — v)/N respectively, one finds that there are 
three types of elementary excitation modes. A hole in 
the A-configuration, /i-configuration, or (^-configuration 
(A-hole, /i-hole, and i^-hole respectively hereafter) creates 
a SU{A) multiplet labeled by the highest weight (1, 0, 0), 
(0, 1,0) or (0,0, 1) respectively. The ^-hole has a six di- 
mensional representation, the A- and i/-holes have four 
dimensional representations. We call these excitations 
the SU{A) flavorons. The energies of these elementary 
excitation are 



£.(A) 
er{v) 



Jtt/2 



V2cosh(A7r/2)-l' 

j7r/2 
cosh (^7r/2) ' 
j7r/2 
V2cosh(p7r/2)-f l' 



(15) 



where A, /i and v stand for the positions of holes in the 
corresponding rapidity configurations. These excitation 



energies vanish when the positions of holes go to infinity 
in the thermodynamic Umit. Therefore they are gapless 
modes. The momenta of the excitations are given by 



Per (A) — 2 tan 

Pu){p.) = 2 tan" 
p^(p) = 2 tan ' 



(\/2+l)tanh(A7r/4) 
tanh ifln/4)] - n/2, 
{V2~ l)tanh(P7r/4) 



3^4, 



7r/4. (16) 



Ehminating the rapidities in eqs. (Eq) and (Jig), we have 

Jtt r 



£a{Pa 
SojIPu. 



2 
Jtt 

'Y 

Jtt 



\/2cos(p^ + 37r/4) + l 
cos{p^ +7r/2), 



^-t{Pt) = -z- V2cos(pr + tt/A) - 1 



(17) 



where p,, e [-3tt/2,0], p^ £ [-tt,0], pr G [-tt/2,0]. 

We are now in the position to relate the elementary 
excitations to the spin and orbital in the original model, 
eq. (0). The quadruplets (1, 0, 0) or (0, 0, 1) are flavorons 
carrying both spin 1/2 and orbital 1/2 with energies e^ 
or St, the hexaplet (0, 1,0) describes flavoron s carrying 
either spin 1 or orbital 1 with energy e^. The spectra of 
these three types of excitations are plotted in fig.(|l]). In 
comparison to Sutherland's results, eq.([l7|) differs only 
in that each mode is shifted by a different constant in 
momentum arising from the tt term in Eq. (S). It seems 
more convenient to use the present version, eq. (|l7|), to 
study the spectra of the magnon-type excitations in the 
subsection V.C. 



(^i{q) = -exp(igAo - \q\/2). 

Our calculation shows a complete cancellation in eq. ( [l9| ) . 
Therefore a 2-string in the A-configuration does not 
change the energy. We also find AJ — N J (T{X)dX + 2, 
where a{X) is the density of the real roots. 

Similar results are found for the 2-strings in the v- and 
/i-configurations. The equations for a z^-configuration are 
given by, 

uJhin) = -Klip - vo)/N, 

Th{v) = [K3/2{V - i^o) + Ki/2{V - V(,)]/N 

and M" = N J T{v)di' + 2. By solving cri(A) = 0, we 
obtain 

l\E = -2ttJ I Ki/2<^i{X)dX = 0. 

And the equations for a 2-string in the /i-configuration 
are given by, 

ah{X)=-Ki{X-ixo)/N, 

uJhifJ.) = [Ki/2{fJ. - A^o) + i^3/2(/^ - l^o)]/N, 

and M' — N J u){ii)d^+2, we obtain iTi(A) — and hence 
A£' — 0. Although these three types of the 2-strings do 
not contribute to the energy, they do contribute to the 
quantum numbers of spin and orbital, and to the highest 
weight of the SU{4) representations. 



B. The complex roots 

Because of the existance of the complex roots p5| in 
the solution set of the Bethe ansatz equations (|7|), we 
must consider their contributions, particularly from the 
2-strings. In this case we need to rederive the integral 
equations (11). We obtain the same equation formally 
but now the inhomogeneous terms ah{X), Uhip), and 
Th{v) include also the contributions from the complex 
roots. A 2-string in the A-configuration, A± = Aq ± i/2, 
introduces additional terms in eq. (pT|). As a result, we 
have 

a,. (A) = [i^3/2(A - Ao) + i^i/2(A - Ao)]/iV, 

Th{u) = 0. (18) 

The energy is given by 

/oo 
Kil2{X)ai{X)dX 
-OO 

-27rJ[/^i/2(Ao + ^/2) + i^i/2(Ao - i/2)], (19) 
and the integral equation leads to 



C. Generalized magnon type excitations 

The flavorons discussed in the previous subsection are 
elementary excitations of the system. These flavorons 
may combine to form composite excitations similar to 
the magnon excitations in the 1-dimensional spin chain 
^6| , which are of interest to experiments and numerical 
simulations. In such a construction, the structure of the 
decomposition of the direct product of the SU{A) funda- 
mental representation must be taken into account. Let 
us consider N = 4n, the decomposition brings about a 
direct sum of a series of irreducible representations, i.e., 
(0, 0, 0), (1, 0, 1), (0, 2, 0), (2, 1, 0), (4, 0, 0) etc. The com- 
posite excitation states include both the singlet (0, 0, 0) 
and the multiplets of 15-fold (1, 0, 1), of 20-fold (0, 2, 0), 
and of 45-fold (2,1,0) or of 35-fold (4,0,0) etc. Those 
multifold excitations are the generalization of magnon 
excitations to the spin systems with orbital degeneracy. 
One A-hole and one z^-hole together create a 15-fold mul- 
tiplet with excitation energy and momentum, 

A£;(i5) =£^(A) + £^(p), 

AP(15) =P<t(A) +Pr{v), 



which is a pair of flavorons of cr-type and r-type. Two 
/i-holes create a 20-fold multiplet with 

AP(20) ^PujiPl) +Plo{P.2)- 

The 45-fold multiplet is a three-hole state created by two 
A-holes and one /i-hole, for which the excitation energy 
and momentum are 

A£^(45) = £^(Ai) + ea{X2) + £uj{p.), 
AP(45) = Pa{Xl) + Pa{X2) + Puj(fi)- 

Four A-holes create a 35-fold multiplet with 



The w-type excitation mode is a moving hexaplet with 



AE, 



(35) 



AP< 



(35) 



4 

E 

4 

= E 



e^lAj), 



P<t(Aj 



The singlet excitation is obtained by placing a A-hole, a 
i/-hole and three 2-strings in A, n and i^-configurations 
respectively. The singlet is degenerate with the 15-fold 
multiplet in energy, i.e. 

In the above equations, Co-IA), eu{p,) and Criv) are given 
by eq.(p^, and Pa(X), pui{p), and Pr(i') are given by 

The energy-momentum dispersion of various magnon- 
types of excitations are plotted in Fig. (0) . In the spec- 
trum calculations, we have used the periodicity in mo- 
mentum P, so that AE{P + 27r) = AE{P). For in- 
stance, for the 45-fold degenerate states, A_E(45)(P) = 
Co-fe) + Ecrfe) + ecuiqs), where P ^ qi + q2 + q^,, with 
modula 2tt. In a recent paper in Ref. 15, the lower-lying 
excitations of model (H) were calculated numerically for 
finite systems. Their results are consistent with ours in 
Fig. 2. In particular, both the numerical calculations 
and the present Bethe ansatz solutions show the following 
feature: as the momentum \p\ increases from to tt, the 
lowest excitations are changed from the 15- fold degener- 
ate states to the 45-fold degenerate states at |P| = tt/2. 

In Fig. 0, we illustrate these generalized magnon 
types of composite low lying excited states. We start 
with a typical configuration of the the ground state in 
Fig. 1^ (a), and the various generalized magnon exci- 
tations are created from the ground state by introduc- 
ing two or more flavorons as shown in (b)-(f), which 
arise from various possible flips of spin or orbital or 
both. The flavoron indicated by the dashed box in 
fig.(0) moves in the background of the S'[/(4) singlet 
carrying both energy and the quantum numbers. The 
CT-type (flavoron) excitation mode is a moving quadu- 

plet with I }_ > being the local highest weight state. 



T i 



> — I ^ J^ > being the local highest weight 



state. The r-type excitation mode is a moving quadu- 
' ^ ^-It£_[> being the 



plet with I 1 i 



i T 



T 



> 



T 



> - 



local highest weight state. 



VI. SUMMARY AND ACKNOWLEDGMENT 



In this paper we have used Bethe ansatz method to 
discuss extensively the ground state and various types of 
the low lying excited states of a Heisenberg spin chain 
with two-fold orbital degeneracy in the limit of S'C/(4) 
symmetry. There are three types of elementary excita- 
tions in the present model. Two of them carry spin 1/2 
and orbital 1/2, and both arc four-fold degenerate. The 
third one carries either spin 1 or orbital 1, and is six- 
fold degenerate. We have constructed magnon types of 
composite excitations and calculated their spectra. 

We thank G.S. Tian for a reference. YQL acknowl- 
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FIG. 1. The spectra of the three type quasi-particles 
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FIG. 2. Th e dispersions for various types of the generalized magnons. E and P stand for the the corresponding /S.E and 
AP in Seclvcl 
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FIG. 3. (a). A < SjTj > staggerily ordered state to demonstrate the ground state, the SU(4) singlet, (b). flavorons of 
one (7-type and one T-type compounded symmetrically to form a 15- fold multiplet, and (c). compounded anti-symmetrically 
to form a singlet, (d). The 20-fold multiplet being compounded of two flavorons of (j;-type. (e). The 45-fold multiplet being 
compounded of two flavorons of cr-type and one flavorons of tj-type. (f). The 35- fold multiplet being compounded of four 
flavorons of cr-type. 



